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ABSTRACT 

The redshifted 21cm line of neutral hydrogen (Hi), potentially observable at low radio frequencies 
(~ 50 — 200MHz), is a promising probe of the physical conditions of the inter-galactic medium during 
Cosmic Dawn and the Epoch of Reionisation (EoR). The sky-averaged Hi signal is expected to be ex¬ 
tremely weak (^100mK) in comparison to the Galactic foreground emission 10"*’K). Moreover, the 
sky-averaged spectra measured by ground-based instruments are affected by chromatic propagation 
effects (^tens of kelvin) originating in the ionosphere. We analyse data collected with the upgraded 
BIGHORNS system deployed at the Murchison Radio-astronomy Observatory to assess the signifi¬ 
cance of ionospheric effects on the detection of the global EoR signal. The ionospheric effects identified 
in these data are, particularly during nighttime, dominated by absorption and emission. We measure 
some properties of the ionosphere, such as the electron temperature (Tg Ri470K at nighttime), mag¬ 
nitude, and variability of optical depth (tiqqmhz ~0.01 and 6t ~0.005 at nighttime). According to 
the results of a statistical test applied on a large data sample, very long integrations (^100 h collected 
over approximately two months) lead to increased signal to noise even in the presence of ionospheric 
variability. This is further supported by the structure of the power spectrum of the sky temperature 
fluctuations, which has flicker noise characteristics at frequencies > 10“® Hz, but becomes flat below 
ft! 10“® Hz. Hence, we conclude that the stochastic error introduced by the chromatic ionospheric 
effects tends to zero in an average. Therefore, the ionospheric effects and fluctuations are not fun¬ 
damental impediments preventing ground-based instruments from integrating down to the precision 
required by global EoR experiments, provided that the ionospheric contribution is properly accounted 
for in the data analysis. 

Keywords: cosmology: observations - dark ages, reionization, first stars - atmospheric effects - meth¬ 
ods: observational - methods: data analysis - instrumentation: miscellaneous 


1. INTRODUCTION 

Detection of the highly redshifted 21 cm signal due 
to the spin-flip transition in neutral hydrogen can be a 
powerful probe of the early Universe {z > 15) and the 
Epoch of Reionisation (EoR; 15 > z > 6) in particu¬ 
lar. Observational efforts are being undertaken on sev¬ 
eral fronts ranging from large interferometric arrays to 
single antenna experiments. The latter aim to detect 
the signature of the EoR in the integrated spectrum of 
the entire visible sky at low frequencies - the “global” 
EoR signal. This measurement is extremely challeng¬ 
ing, primarily because the expected signal (~ 100 mK) is 
4 — 5 orders of magnitude weaker than the foregrounds 
(~1000K). The global EoR approach has already re¬ 
sulted in a lower limit (Az > 0.0 6) on the duration of 
the EoR (|Bowman fc Roger^l2010t) obtained by the Ex¬ 
periment to Detect the Global EoR Step (EDGES). 

Achieving precision on the order of a few milliKelvin 
requires an extremely stable radiometer system and a 
nearly perfect understanding and calibration of the sig¬ 
nal. For a given frequency, the required integration time 
depends on the receiver noise temperature and frequency 
bin width selected to optimize data analysis (limited by 
the resolution of a particular spectrometer). In princi¬ 
ple, according to estimates for the Broadband Instrument 


for Global Hydrogen Reionisation Signal (BIGHORNS; 
Sokolowski et al. 2015, hereafter referred to as S15) sys¬ 
tem with a receiver noise temperature ^90 K, a statisti¬ 
cal error of ^10 mK in a 1 MHz bin at frequency 80 MHz 
can be achieved in ^4 hours of integration (^ImK re¬ 
quires ~400 hours) at optimal, “cold” sky, conditions (i.e. 
nighttime when the Galactic Genter is below the horizon 
and sky noise temperature at 80 MHz is ftl200K ). In 
practice however, due to the uncertainty of the under¬ 
lying sky spectrum, which needs to be fitted and sub¬ 
tracted, integration times longer by as much as an or¬ 
der of magnitude are required to achieve residuals of 
~10mK. Moreover, precision of ~10mK might not be 
sufficient for data analysis methods to constrain reioniza¬ 
tion models and recover physical parameters (represented 
by predicted turning points in the sky-averaged signal) of 
the tr ansition from the Dark Ages to the re ionized Uni¬ 
verse (jBiirns et al.ll201^ iHarker et al.ll2012[l . Therefore, 
integrating down to ~mK precisions requires averages 
over data collected during multiple nights. 

This basic estimate, however, does not account for 
additional stochastic erroiQ which may come from in- 

^ In order to reflect the nature of ionospheric fluctuations, we will 
further call the additional error which they introduce “stochastic 
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strumental or other effects such as propagation of the 
sky signal through the Earth’s atmosphere. The atmo¬ 
sphere can affect the cosmic signal propagation via sev¬ 
eral effects. The most important for the global EoR sig¬ 
nal experiment are chromatic effects in the ionosphere. 
Achromatic refraction in the troposphere is less impor¬ 
tant because it does not affect the frequency structure 
of the signal. Moreover, its variability (due to varia¬ 
tions in ambient temperature and atmospheric pressure), 
which we estimate to be <2 K above 50 MHz, is at least 
one order of magnitude smaller than variability of the 
chromatic ionospheric effects. Nevertheless, the tropo¬ 
spheric refractio n may have an indi rect impact through 
duct formation (|Hitnev et al.lfTOSSf) which enables long 
distance radio-frequency interference (RFI) propagation 
and therefore occasionally causes significant data exci¬ 
sion rates. In the following subsections we will present 
some properties of the ionosphere most relevant to the 
global EoR experiments. 


1.1. Ionosphere 

The ionosphere is the part of the Earth’s atmosphere 
where there is sufficient ionization to affect the propaga¬ 
tion of radio waves. It consists of several layers: D-layer 
at height 50 — 90 km; E-la yer at 90— 140 km; and E-layer 
at he ights above 140 km (|Davieslll990t lEvans fc Hagford 
1196811 . The most relevant to single antenna, high preci¬ 
sion measurements of sky brightness temperature at low 
frequencies are the D-layer and F-layer. 

(i) F-layer: most of the ionospheric electron column 
density is accounted for by the F-layer, where electron 
density is ~ 10^^e“/m^ during the daytime and typ- 
ically an order of magnitude less du ring the nighttime 
(jDavie^ll990l: lEvans fc HagforsI 1 19681 1. Measurements of 
the total electron content (TECj^l, routinely monitored 
by measuring the propagation delay of GPS signals, are 
dominated by electrons in the F-layer. 

Due to very high electron content the F-layer domi¬ 
nates refractive effects, which may play a non-negligible 
role for global EoR experiments by shifting radio sources 
within the antenna beam. Based on equation 6 in [Bailed 
(fl^ . we estimated the average position offsets of radio 
sources for typical nighttime electron density 10^^ e~ jrr? 
to be ^ 0.1 arcmin at 100 MHz and zenith angle Q = 45° 
(corresponding to ^ 0.3 arcmin at 1 00 MH z and zenith 
angle 60° cited bv iThomnson et HI (|2007tl l. According 
to the Bailey’s equation, the expected intra-day vari¬ 
ations of the average offset are proportional to varia¬ 
tions in_dectron density of typically ^20-30% (see Sec- 
tion 13.5.21). Using the data from the GLEAM survey 
(Wavth et al.l j201,^ performed with the MWA telescope 
(|Tingav et al.l 1201^ located at the Murchison Ra dio- 
astronomy Observatory (MRO), lArora et al.l (|2015[ 1 re¬ 
port nighttime offsets of the source positions of the or¬ 
der of a fraction of arc min at 100 MHz and scaling as 
1/z/^. On the other hand. lLoi et al.l(|2015f) report slightly 
higher ^1 arcmin source offsets at 100 MHz due to elec¬ 
tron density gradients (plasma tubes) in the ionosphere, 
but such prominent effects were only observed occasion- 


error” . 

^ TEC is the total number of electrons per rnl (electron col¬ 
umn density) along a given line-of-sight (typically determined for 
vertical direction). 1 TEC unit is a column density of 10^® e~I’m? 


ally and therefore can be treated as a non-typical behav¬ 
ior. 

On the other hand, signal loss in the ionosphere 
measured in dB is expected to be of the order of 
Ldn ~ iJcTECLiWO/i^MHz)'^ (jEvans fc HagforsI 

119681) . where i7c is the mean rate of electron collisions, 
and TECl is the mean total electron content in the layer 
L, expressed as a product of mean electron density (iVe) 
and the width of the ionospheric layer (Aft.^). Therefore, 
due to low collision rates ^ 10^ Hz (|Thrane fc PiggottI 
119661: lAggarwal et al.lll97^ . nighttime absorption in the 
F-layer is considered to be smaller than in the D-layer; 
but due to high e lectron density it can be s ignificant dur¬ 
ing the daytime (jSteiger fc Warwi^ll96lD . 

(ii) D-layer: is the lowest part of the ionosphere 

where electron densities are ^2-3 orders of magnitude 
lower than in the F-layer. High electron densities in 
the D-layer exist only during the daytime due to in¬ 
solation. Typical daytime electron density in the D- 
layer is ~ 10®e“/m^, falling to a residual density of 
^ 10®e“/m^ at night. Due to significant collision rates 
of the order of Vr. ^ 10^ Hz (iVuthaluru et al.l 120021 : 
lAggarwal et al.l 119791 : iDaviesI I1990D the D-laver dom¬ 
inates nighttime ionospheric absorption. The night¬ 
time residual electron density causes absorption ef¬ 
fects of order Lds ^ ^ {IQQIvmhz Y _A 

0.01 dBx (ICj O at zenith (lEvans fc: Hagfoielll968l : 
IDaviesI [19901 : IThomnson et al.ll2007D . which may have a 
significant impact on global EoR measurements. The 
magnitude of typical nighttime absorption at 100 MHz, 

~ 0.01 dB, corresponds to optical depth of the 
ionosphere r Ri 0.0023. 

(iii) Thermal em ission: the ionosphere emits ther¬ 
mal radio emission (jHsiehl I1966D with brightness tem¬ 
perature Tiono ^ TTe, where r is the optical depth of 
the ionosphere and is the electron temperature. Be¬ 
cause Tiono is proportional to optical depth it is expected 
that nighttime emission is dominated by the D-layer, 
whilst during the daytime there should be a significant 
contribution from the F-layer. The electron tempera¬ 
ture in the ionosphere has been measured with a variety 
of me thods ranging from rocket measurements (jOvamal 
120061 ) to ground-based observations. The electron tem¬ 
peratures of the lower ionosphere (D-layer) measured at 
frequencies 1.6 and 2 MHz, at which radio signals from 
outer space are entir ely reflected, are ty pically in the 
range ^ 240 — 290K ([Pawsev et al.l[T95ll) . The electron 
temperature s at hi gher altitudes (F-layer), measured by 
iZhang et al.l ()2004[ ) using incoherent scatter radar (ISR) 
experiments, are typically ^ 800 — 1000 K at 125 < h < 
175 km and ~ 1000 — 2000 K at 275 < h < 325 km, and 
exhibit seasonal and solar cycle variations. 

Due to these effects the ionosphere, besides affecting 
ground based and satellite radio-communications, signif¬ 
icantly influences low frequency radio signals from outer 
space observed on the ground. Moreover, the properties 
of the ionosphere vary in time on multiple timescales due 
to factors affecting the amount of ionization such as, for 
example, the solar cycled, solar activity on a particular 
day, or cosmic-ray flux. All these effects have resulted 

® Beginning of 2015 is just after the maximum of cycle 24 
(http://solarscience.msfc.nasa.gov/predict.shtml) 
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in an on-going discussion in t he literature on the im- 
portance of ionospheric effects (|Vedantham et aLll2014f) 
and whether the global E oR signature can feasibly be 
detected from the ground (iDatta et al.ll201!Tll . 

1.2. Ionospheric impact on the global EoR detection 

The impact of the aforementioned ionospheric effects 
on the detecti on of the global E o R has rece ntly been in¬ 
vestig ated by iVedantham et al.l (|2014f) and iDatta et al.l 
(|2015ll . These studies have shown that ionospheric ab¬ 
sorption and refraction can significa ntly affect the fre - 
quency structure of the cosmic signal. IDatta et al.l (120151 1 
conclude that long timescale variability of ionospheric 
absorption and refraction resulting from observed TEC 
variability may impede even perfect ground-based instru¬ 
ments from integrating down to the required ~mK pre¬ 
cision. Hence, it is vital to identify ionospheric effects 
in the sky-averaged radio spectrum at the frequencies 
of interest (50 — 200 MHz) and assess their significance 
for global EoR experiments. The observation of iono¬ 
spheric effects at frequencies 80 — 185 MHz in the sky- 
averaged radio spect rum have recently b een reported by 
the EDGES group (|Rogers et al.ll201^ . The EDGES 
experiment measured variability of absorption (optical 
depth) and magnitude of emission (electron temperature) 
of the ionosphere, which proves that these effects do play 
an important role and must be accounted for in data 
analysis. 

Additionally, such measurements present an interest¬ 
ing alternative technique to measure the properties of 
the ionosphere at low radio frequencies with wide-band 
total power radiometers. Typically, relative absorption 
is measured by narrow-band relative ionospheric opacity 
meters (riometersfl at l ower frequencie s of ^ 30 MHz 
(iLittle fc LeinbacH ~ fl959ll or 18.3 MHz (|Mitra fc ShainI 
1195311 . Ha.s'^ on multi-day observations, a so-called 
‘quiet day curve’ (QDG) of the maximum power ob¬ 
served by the riometer during a state of very small, 
residual ionos pheric absorption (typical l y dur ing the pre¬ 
dawn hours (IBelikovich fc BenediktovHl969ll l is deter¬ 
mined and re lative absorption is measu red with respect 
to the QDG (|Krishnaswamv et al.lll98^ . Similar tech¬ 
niques can be applied to wide-band sky-averaged radio 
spectra from global EoR experiments, which makes them 
useful tools for ionospheric studies. Moreover, wide-band 
data from these instruments enables analyses of the fre¬ 
quency structure of ionospheric absorption and emission, 
which are not available from traditional narrow band ri- 
ometers. Finally, applications of the riometer approach 
are not limited to ionospheric studies, as they turned 
out to be useful tools even for inferring X-ray flux of the 
large st solar X-ray flare ever recorded (|Brodrick et al.l 
[ 20 (^ . 

This paper presents an analysis aiming to assess the 
significance of ionospheric effects on the detection of the 
global EoR signal. Data were collected in 2014-2015 

^ Throughout this manuscript we use the term “optical depth” 
for the measured absorption. However, we find the naming in the 
riometer field confusing as people typically talk about “opacity” 
and even the term riometer is an acronym for a relative ionospheric 
“opacity” meter, whilst in fact everybody seems to measure what 
is strictly speaking “optical depth”. Probably, like in many other 
fields this is some historical relic, but we did not get to the bottom 
of the problem. 


by the BIGHORNS system at the Murchison Radio- 
astronomy Observatory (MRO) in Western Australia. In 
Section [2] we describe the BIGHORNS system as it was 
deployed at the MRO in October 2014, mainly focusing 
on modifications with respect to the system presented in 
S15. In Section [3] we present data processing and data 
analysis leading to identification of the ionospheric effects 
in the data. We compare these results to the expected 
contributions from ionospheric absorption/emission and 
refraction. We also present some properties of the iono¬ 
sphere, such as electron temperature, relative optical 
depth, and its variability derived from the data analysis. 
In Section [5] we compare our derived electron tempera¬ 
ture with model predictions. Finally, in Section[5]we per¬ 
form a statistical test to verify whether the observed data 
can be integrated down to the required ~ mK precision 
by averaging large numbers of spectra in the presence of 
a variable ionosphere. We compare observed variability 
in the real data with the variability expected for diurnal 
changes in antenna temperature due to Earth rotation. 
In order to better understa nd the result s of th is test in 
the context of the results of IDatta et al.l ()2015r i. we also 
study the power spectrum of observed fluctuations in the 
sky temperature. 

2. INSTRUMENT DESCRIPTION 

The detailed description of the initial installation of 
the BIGHORNS total power radiometer can be found 
in S15. The BIGHORNS system in a very similar form 
was deployed at the MRO in October 2014. The major 
difference with respect to the system described in S15 
is the antenna. The portable biconical antenna, which 
was used during short remote deployments, was replaced 
with a conical log-spiral (GLS) antenna (Fig. [T]). The 
antenna was painted white in order to protect the cop¬ 
per conductor from oxidizing and to protect it from bad 
weather conditions. The balun is installed at the apex 
of the cone under a protective cap and is connected to 
the front-end receiver via a cable of approximately 2.5 m 
length. 

The design details of the antenna and m otivation for 
the upgrade can be found in Lliwanil (120131 1 and S15. In 
summary, the upgrade was motived by the following rea¬ 
sons. Firstly, the response (beam pattern and the re¬ 
flection coefficient) of the GLS is, by design, more fre¬ 
quency independent than the response of the previously 
used biconical antenna. Secondly, it is very well matched 
(<-15dB) across a wide frequency band (~70-180MHz), 
which improves signal to noise ratio and facilitates a more 
accurate calibration over a wider frequency band. Par¬ 
ticularly, the GLS antenna significantly improves sensi¬ 
tivity at low frequencies (70MHz "^v 150 MHz), which 

are of main interest for the EoR experiment. Neverthe¬ 
less, its reflection coefficient has more frequency struc¬ 
ture (i.e. ripples or spikes) compared to the biconical 
antenna (Fig. [2]). This frequency structure, which re¬ 
sults from the constructive and destructive interference 
of the RF-signals reflected from the unterminated ends of 
the spiral arms, may pose a problem for the data analysis 
looking for a tiny (~ mK) global EoR “ripple” amongst 
the foregrounds (assumed to be smooth as a function 
of frequency). The comparison of the reflection coef¬ 
ficients of the biconical and GLS antennas, which were 
measured at the end of the cables connecting the balun to 
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Figure 1. The conical log-spiral antenna as deployed at the MRO 
in October 2014. The antenna was painted in order to protect the 
copper conductor from the weather. 


the front-end receiver, is presented in Figure O Thirdly, 
the exact knowledge and accessibility of all its compo¬ 
nents (such as the balun, output cable etc.) enabled 
us to better characterize the CLS antenna. Finally, the 
CLS antenna is a left-hand circularly polarized (LHCP) 
antenna which slightly helped in suppressing reception of 
undesired right-hand circularly polarized (RHCP) signals 
from the ORBCOMM satellites. We estimate that the 
CLS antenna attenuates the undesired ORBCOMM sig¬ 
nals (at 137-138 MHz) by a few dB (see also Section [SH]). 

The antenna was modeled in the FEKO simulation 
package. Its antenna pattern is within ~10% frequency 
independent and symmetrical in azimuthal angle (</>) in 
75-200MHz band (Fig. [S]). Hence, we will further ap¬ 
proximate its response as frequency and azimuthal angle 
independent. Upgrading to the better-matched antenna 
allowed us to remove the 3 dB attenuator which was orig¬ 
inally installed between the biconical antenna and the 
front-end box to suppress the effects of the impedance 
mismatch between the two. Therefore, the receiver noise 
temperature of the new system (without any additional 
attenuation between the antenna and front-end) was re¬ 
duced to Trcv ~180K. 

3. DATA ANALYSIS 

The system collected data in exactly the same mode as 
described in S15. Therefore, the data acquisition param¬ 
eters will be only briefly reviewed here. In the present 
configuration the system uses cold calibration schema 
and switches between the antenna (for fa 15.5 sec) and 
a reference source (for 5 sec), which is a 50 H termina¬ 
tor at ambient temperature (measured with a portable 
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Figure 2. Comparison of the reflection coefficient magnitudes 
(upper figure) and phases (lower figure) of the log spiral antenna 
and previously used biconical antenna. 


temperature logger). The data were collected continu- 
ously with 100% duty cycle (no gaps in tim(|3) at ap¬ 
proximately 50 ms time and 117.2 kHz frequency resolu¬ 
tion. The FPGA based spectrometer measured the spec¬ 
tra in 0 — 480 MHz band and the analog filters limited 
the bandwidth to ~ 40 — 350 MHz. The power spectra 
produced by the s pectrometer were saved to FITS hies 
(iPence et al.ll2010tl , which consisted of interleaved blocks 
of antenna and reference data. Each FITS hie consisted 
of 1000 integrations, thus every hie covered a 50 sec in¬ 
terval. 

3.1. Initial data pre-processing and calibration 

The details of initial data pre-processing can be found 
in S15. In summary, it consisted of the following main 
steps: 

1. State identification - The calibration RF-switch 
operates autonomously. Therefore, the first data 
processing step identifies integrations collected on 
the antenna and reference source and saves the lo¬ 
cation of the transitions to FITS file headers. 

2. RFI excision - The entire integration was ex¬ 
cised if the power measured by the spectrometer ex¬ 
ceeded a threshold of -35 dBm (ai 3 x 10“"^ mW)0 in 
any single channel or the sum of power in the entire 
band was above the median galactic noise level plus 
the threshold of-37.5dBm (si 1.8xl0-4mW). The 
data accepted by these two criteria were flagge d 
by the AO-flagger software (jOffrinea et al.l[201^ . 
in order to identify individual channels affected by 
RFI. Finally, the data were checked with variabil¬ 
ity criteria (see Sec. 13.211 in order to exclude sig- 


® The recorded data are continuous, but the sky data (collected 
on the antenna) have 5 sec gaps due to the aforementioned switch¬ 
ing to a reference source. 

® dBm is the ratio of measured power referenced to 1 mW and 
expressed in decibels, PdBm = 10 ^ogio(Pmw) 
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Figure 3. The antenna pattern of the conical log spiral antenna 
simulated in the FEKO software. The upper figure shows the 
simulated CLS antenna gain as a function of elevation angle at 
at 90, 110, 150, and 190 MHz (red lines) with a Gaussian beam 
gain~exp(—(6/52°)^) (blue line). The lower figure shows antenna 
gain as a function of azimuthal angle at elevation angle 70°. The 
antenna response is within ~ 10 % frequency and azimuthal angle 
independent. 

nificantly variable data (due to solar activity or 
lightning). The average excision rate due to high 
power (>-35 dBm) from the ORBCOMM satellites 
(at 137-138 MHz) of the system with the CLS an¬ 
tenna (without the 3 dB attenuator between the 
antenna and the front-end) was fvO.7% comparing 
to ~10.9% with the previously used biconical an¬ 
tenna (after adding 3 dB due to an additional atten¬ 
uator). Therefore, based on shift of the two power 
distributions, we estimate that the CLS antenna 
suppresses the undesired ORBCOMM signals by a 
few dB (^3-5dB). 

3. Data reduction - The data unflagged by the RFl 
excision procedure were reduced to a more manage¬ 
able volume by averaging (with a median estima¬ 


tor) every Nant = 700 integrations on the antenna 
accepted by the RFl excision criteria and ACe/ in 
the range 200 — 300 on the reference load. After 
this step the data formed an averaged dynamic 
spectrum where each integration was an average 
of Tant = 700 X 0.05 = 35 sec of single integrations 
in the case of antenna and Tref ~ 12.5 sec in the 
case of reference source. The dynamic spectra from 
the antenna and reference were saved into separate 
FITS files. To enable correct error propagation, 
the numbers of averaged individual 50 ms integra¬ 
tions were saved to FITS files (also separate for 
the antenna and reference data) for every averaged 
integration and every frequency channel. Hence, 
calculation of the statistical error on the averaged 
spectra used actual integrations times (after exci¬ 
sion of RFl affected data). The resulting reduced 
data were split by calendar months and saved into 
separate FITS hies representing dynamic spectra 
for each calendar month. Further data analysis was 
performed on the reduced dataset. 

4. Calibration - each averaged antenna integration 
in the dynamic spectrum of a given month was cali¬ 
brated by the corresponding averaged reference in¬ 
tegration, according to the calibration procedure 
described in S15. The procedure uses the power 
measured on the reference source and information 
from the ambient temperature probe to calibrate 
the antenna signal in units of temperature; next, 
the rehection coefficients (of the antenna and the 
front-end receiver) are used in order to correct for 
the reflection of the sky signal from the input of 
the receiver; hnally, the receiver noise temperature 
and the input noise of the receiver reflected back 
and forth between the antenna and the receiver (so 
called reflected noise) are subtracted. 

The only major difference in the calibration pro¬ 
cedure was that the reflection coefficient of the 
conical log spiral antenna measured at the MRO 
(Fig. [2]) was used instead of characteristics of the bi¬ 
conical antenna and there was no additional atten¬ 
uation between the antenna and the front-end. As 
a consistency check, the ca librated data were com¬ 
pared with a sky model (|de Oliveira-Costa et al.l 
12001 integrated with an antenna pattern simu¬ 
lated in the FEKO electromagnetic simulation soft¬ 
ware. The relative difference between the predicted 
antenna temperature from the model and the cal¬ 
ibrated data is approximately 10%, which is con¬ 
sistent with uncertainties in the model tempera¬ 
ture scale. The resulting calibrated dynamic spec¬ 
tra and corresponding statistical errors (calculated 
according to equation 13 in S15 for appropriate in¬ 
tegration times on the antenna and reference) were 
saved to separate FITS files. The contribution of 
the statistical error of averaged reference integra¬ 
tions (jref ~ 12.5) to the overall budget of the 
statistical error is ~50%, and may be improved 
in the future in order to reduce the required in¬ 
tegration time. Therefore, after calibration, two 
fits files containing the dynamic spectrum and its 
error were created for each analyzed month (Octo- 
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ber/November 20143, December 2014, and January 
2015). 

Further data analysis presented here was performed 
separately on each calibrated monthly dynamic spectra. 


3.2. Excision of significantly variable data 

The data processing described in S15 solely focused 
on nighttime data without any consideration of daytime 
data, which is regarded as low value for the global EoR 
analysis. However, in the presented analysis we aim to 
derive some properties of the ionosphere as a function 
of local time including the daytime. Therefore, in order 
to include only quiet-Sun data and exclude data affected 
by solar activity we needed an effective way of excluding 
data affected by variations on timescales of seconds and 
less. Some part of these data are rejected by the AO- 
flagger software in the RFI excision process, but in order 
to further clean the data sample we developed a pro¬ 
cedure which identifies and flags data affected by short 
timescale variability. The procedure acts on individual 
FITS files (1000 integration each) at 50 ms resolution af¬ 
ter the state identification procedure and is based on the 
radiometer equation: 


^expiy') 


PH ^ 

"x/ P^int 


( 1 ) 


where PH is the power as a function of frequency 
measured by the radiometer system in arbitrary units, 
CexpH is the standard deviation of the measured power, 
B is the resolution of the frequency bin, and Tint is the 
integration time. In the absence of short timescale vari¬ 
ability, such as solar activity or lightning, the expected 
standard deviation (aexp) as a function of frequency can 
be calculated from the mean or median spectrum accord¬ 
ing to equation [H On the other hand, the actual stan¬ 
dard deviation (aobs) can be calculated from a sample of 
antenna integrations. 

In the presented analysis we calculated the observed 
standard deviation from all the antenna integrations in 
a single FITS file (typically ^ 700x50ms integrations). 
The main goal of this criterion was to identify and excise 
any kind of short timescale variability (< 35 sec) such as 
solar activity (second and even millisecond timescales) 
or lightning (even sub-millisecond timescales). Longer 
timescale variability (^35 s) was most likely not identi¬ 
fied by this procedure. Nevertheless, lack of significant 
variability on longer timescales was verified on the dy¬ 
namic spectrum normalized by a median 24 h (such as 
Fig. 18 in S15). The same integrations were used to 
calculate the mean power spectrum PH ^-nd the ex¬ 
pected standard deviation aexpH calculated accord¬ 
ing to equation [TJ Under nighttime conditions and in the 
absence of any short timescale variability (e.g. nearby 
lightning) the expected and observed standard deviation 
curves lie on top of each other (Fig.|4]). However, in the 
presence of solar flares or thunderstorms the observed 
standard deviation is signihcantly larger than calculated 
according to equation [T] (Fig. |4]). Therefore, the differ¬ 
ence between the observed and expected standard devi- 


^ In October 2014 only a few days worth of data were collected. 
Therefore, these data were analyzed together with November 2014 
data (forming a 38day dataset). 



Figure 4. Example of difference between observed (black and red 
points) and expected (magenta and blue dashed lines) standard 
deviation of 35 sec data sample. The black points represent the 
observed standard deviation of the nighttime data sample collected 
on 2015-01-16 at 04:00:05 AWST and the magenta dashed line lies 
on top of these data over almost entire bandwidth (except RFI af¬ 
fected channels above 240MHz). The red points, representing the 
observed standard deviation of the daytime data sample collected 
on 2015-01-15 at 13:47:44 AWST in the presence of a solar flare, 
lie much above the dashed blue line corresponding to the standard 
deviation expected for these data. At that time the Galaxy was 
transiting during the daytime, therefore the expected standard de¬ 
viation for daytime (blue long-dashed line) is above the standard 
deviation expected for the nighttime (magenta short-dashed line). 


ation is a good metric of any short timescale variability 
in the data (Fig.[5|). 

The maximum number of consecutive channels Nmax 
(corresponding to a maximum affected bandwidth Bmax) 
in a frequency band 60—300 MHz where the difference ex¬ 
ceeded a frequency independent threshold Lyar = 5 x 10^ 
(corresponding to 5cr around frequency 85 MHz) was used 
as a metric to identify broadband short timescale vari¬ 
ability. The threshold was frequency independent in or¬ 
der to simplify the procedure. However, the method can 
be further refined and made more sensitive by using a fre¬ 
quency dependent value of the Lyar threshold. Several 
frequency bands, for example 88-108MHz (FM band), 
133.9-134.4MHz (aeronautical band), 137-138MHz (OR- 
BCOMM satellite communication band), 242-272MHz 
(satellites), and multiple other narrow (^IMHz) bands 
known for RFI were excluded from the procedure. The 
maximum frequency bandwidth, Bmax , affected by vari¬ 
ability exceeding the Lyar threshold as a function of lo¬ 
cal time is shown in Figure [S] In the later data anal¬ 
ysis, the data were considered to be affected by short 
timescale variability if Nmax for given time exceeded 
a threshold of 200, which corresponds to a bandwidth 
of Bmax ~ 24 MHz. We have tested several values of 
the algorithm parameters and found that Nmax = 200 
{Bmax ~ 24 MHz) is optimal to suppress false identifi¬ 
cations (false positives) due to noise fluctuations in con¬ 
junction with RFI and retain a reasonable amount of 
daytime data (day dependent). The broadband require¬ 
ment (at least Nmax channels affected) together with the 
threshold Lyar > 5cr ensure that the probability of false 
positives due to noise fluctuations is nearly zero. Never¬ 
theless, in order to improve daytime data analysis, the 
procedure can be fine-tuned in order to improve its sen¬ 
sitivity to low-level variability. 

3.3. Data processing 

Our analysis of the ionospheric effects in the sky- 
averaged radio spectrum is based on an approach similar 
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Figure 5. The difference of observed and expected standard de¬ 
viation of 35 sec data sample normalized by the threshold {L^ar)- 
The red dots represent the data collected on 2015-01-15 at 13:47:44 
AWST in the presence of a solar flare and the black triangles repre¬ 
sent the nighttime data collected on 2015-01-16 at 04:00:05 AWST. 
The data points above the blue horizontal line are considered as 
affected by variability. 



Figure 6. Example of variability metrics (maximum affected 
bandwidth) as a function of local time for a sample of data collected 
between 2014-11-01 20:00 and 2014-11-03 20:00 AWST. Clearly the 
nighttime data are very quiet and daytime data show significant 
variability mainly due to solar activity. Occasionally nighttime 
data (not show in this figure) were also significantly affected by 
nearby lightning. The data exceeding the variability bandwidth 
threshold of Ri24 MHz (marked in the image with a red line) were 
excluded from further analysis. 


to riometer observations of relative absorption. In princi¬ 
ple, nighttime spectra collected at a given local sidereal 
time (LST) in the absence of nearby lightning, instru¬ 
mental, or propagation effects should be the same for 
each sidereal day. Therefore, any variation of the same 
LST sky spectra between days is an indication of vari¬ 
ability due to instrumental or local propagation effects. 

The problem becomes more complex for daytime data 
because of the Sun’s contribution. Even the quiet Sun in¬ 
troduces small changes to sky-averaged spectra, because 
its position at a given LST changes slightly every day and 
differs by ^ 30 deg between the first and the last day of 
the month. Thus, for a given LST, the solar position 
within the antenna beam may be significantly different 
and therefore the solar contribution to the sky-averaged 
spectrum is also different. The effect was estimated to 
be of the order a few Kelvins (depending on Sun’s ele¬ 
vation) and was not taken into account in the presented 
analysis, which may have a small impact on daytime data 
analysis. Solar activity events were excised from the an¬ 
alyzed daytime data sample using the methods described 
in Section [221 However, low level solar activity could po¬ 
tentially remain unflagged and contaminate the daytime 
data sample. A more detailed subtraction of the solar 


contribution could further improve analysis of daytime 
data and is planned in the future. 

Calibrated dynamic spectra from each calendar month 
were analyzed separately and integrations accepted by 
variability criteria (Section 13.21) were averaged in I hour 
LST bins (0 — Ih, 1 — 2h, 2 — 3h and so on) forming a 
dynamic spectrum at I h resolution. In order to exclude 
outlier points due to residual RFI, the median estima¬ 
tor was used to average individual Ri35sec integrations 
within 1 h LST bins. Thus, after this step the dynamic 
spectrum consisted of ^ 24 x Ndays individual spectra, 
where Ndays is the number of days in the analyzed cal¬ 
endar month (Fig. |71). 

3.4. Data sample cleaning 

For each 1 h LST bin some additional characteristics, 
such as the mean and standard deviation of the ambient 
temperature and the maximal variability metric {N^ax) 
of the individual integrations within the I h LST bin, 
were also calculated. Based on this information, the fol¬ 
lowing quality criteria were imposed in order to select 
only a consistent sample of I h LST bins not affected by 
instrumental or environmental effects: 

1. Ambient temperature 

The FEKO software simulations indicated that di¬ 
urnal variations in ambient temperature may af¬ 
fect the instrumental response - mainly antenna ef¬ 
ficiency. In the analyzed period, Oct 2014 - Jan 
2015, the ambient temperature measured at the 
MRO varied in the range 285 K - 320 K. There¬ 
fore in order to suppress temperature-related in¬ 
strumental effects to a minimum, only 1 h LST 
bins collected at ambient temperature Tamb within 
Thest ± 5 K were selected for further analysis, where 
Thest is an ambient temperature which maximizes 
the number of spectra in the sample (different for 
every LST time). Typically a few (<7) Ih LST 
bins were excluded for every analyzed month. 

2. Minimum integration time per LST bin 

For each 1 h LST bin in a day’s data, the total in¬ 
tegration time on the antenna Tih was calculated 
as Tih = N X 35 sec, where N is the number of 
35 sec integrations. It is expected to be less than 
1 h because «26% of time was spent on a reference 
source and some fraction of the data were excised 
due to solar activity or RFI. Based on the distribu¬ 
tion of Tih for all 1 h LST bins, we determined that 
Tih = 2500 s with a standard deviation RilOOs. 

The sky temperature changes significantly within a 
Ih interval (as the Galaxy moves across the sky). 
Thus, large gaps in the data were undesired because 
they caused the averaged spectrum of 1 h to be sig¬ 
nificantly different than expected for a fully sam¬ 
pled 1 h period. For example, assume the Galactic 
Center was rising within a given Ih LST bin and 
therefore the sky-averaged signal increases during 
that hour. Further, assume that all the spectra 
in the first half of the LST bin were excised due 
to RFI or variability criteria. Then the resulting 
mean would be biased towards higher sky temper¬ 
atures and therefore not be comparable to other 
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Figure 7. Dynamic spectrum of Oct/Nov 2014 (^SS days) data in Ih LST resolution (upper figure) and a median 24 h day (lower figure). 
The median dynamic spectrum is stretched in the time axis. Each Ih LST bin results from the median of individual 5^35 sec integrations 
accepted by RFI and variability criteria, and magenta (zero value) marks data excised by these criteria (for example in the FM-band around 
100MHz). The 120 MHz channel was entirely excised due to self-generated interference in the spectrometer. The horizontal magenta lines 
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fully sampled bins at the same LST. Hence, in or¬ 
der to exclude 1 h LST bins resulting from averag¬ 
ing spectra with significant gaps, we required that 
Tih > 2200 s, which corresponds to at least 63 in¬ 
dividual 35 sec integrations, collected in any single 
day’s data. Bins with smaller integration time were 
excluded from further analysis. 

3. Sunrise and sunset 

In order to compare only consistent samples of LST 
bins, we excluded bins that contained either a sun¬ 
set or a sunrise as indicated in Table [TJ It is evi¬ 
dent from Table [T] that a given 1 h LST bin could 
be collected during the daytime in one part of the 
month but during the nighttime at another, due 
to the 2h drift of LST at sunrise/sunset during 
a month. We separated these 1 h LST bins into 
a daytime set, and a nighttime set, depending on 
when they occurred. Only one of these sets was 
accepted for further analysis, with the preference 
given to the nighttime set if it contained at least 
four bins deemed acceptable according to other cri¬ 
teria. 

After all the above criteria were applied the number 
of 1 h LST bins was reduced from Ndays to a typically 
smaller number Nacc, where Nacc varied between a few 
(for near sunset/sunrise bins) and 32 (for nighttime bins 
in the Oct/Nov 2014 dataset). Next, for every Ih LST 
spectrum (0—Ih, I — 2h, 2—3h and so on) a reference sky 
spectrum was calculated as a median of the Nacc spectra 
accepted for this LST range (Fig. flower figure). Then 
for each spectrum in the dynamic spectrum at I h res¬ 
olution a corresponding reference median spectrum was 
subtracted. The effect of such a subtraction at the origi¬ 
nal frequency resolution of II7.2kHz for 0 — 1 LST bin in 
November/October is presented in Figure [8] where day- 
to-day variations in received power at the lower frequen¬ 
cies are clearly visible. The differences in the observed 
sky-averaged spectrum in relation to possible ionospheric 
effects will be discussed in the next section. 

Finally, for the sake of fitting ionospheric parameters, 
every 40 channels (~5 MHz) were averaged and the er¬ 
ror of each average channel was calculated as a standard 
deviation of the individual 40 channels. The standard de¬ 
viation of the median-substracted daily signal in several 
5 MHz frequency bins, for all 1 h LST bins, is shown in 
Figure ini We will further refer to this standard deviation 
as (Jiono ■= cr{T^ky ~ Figure [9] demonstrates 

much higher variability in the day-to-day signal at lower 
frequencies over the range of nighttime LSTs, which we 
attribute to ionospheric effects. 

3.5. Expected ionospheric effects 

In this subsection we aim to derive an expected fre¬ 
quency dependence, estimate and compare magnitudes 
of the ionospheric absorption/emission and refraction. In 
the absence of propagation effects due to the ionosphere, 
the antenna temperature, at a given time t, as a func¬ 
tion of frequency resulting from the integration of the 
sky brightness temperature, T^{9,(j)), can be expressed 
by the following formula (equation 1 in S15): 

= P.(0,0)r.(0,<^)dL!, 

M) Jatv 


where P,y{0,(j)) is the antenna pattern and the coordi¬ 
nates {0, (j)) are angles in the antenna frame, which can 
be directly translated to horizon coordinates (azimuth, 
zenith angle) according to orientation of the antenna in 
the field, and Pq = is the normalization 

factor. The above formula is time dependent through 
the implicit time dependence of T^{6,(p). Many equa¬ 
tions and quantities in this paper are time dependent, 
but their explicit time dependence is usually not shown 
for brevity. 

The propagation effects in the ionosphere modify the 
outside-atmosphere sky brightness temperature T^{9,(j)) 
to TP{9,(1)) at the antenna, and therefore the antenna 
temperature observed by a ground-based instrument can 
be expressed as: 

Tant{>y) = ^ [ Pff9,(l))TP{9,(j))dn. (3) 

M) JAtv 

Tant {v) i s obtained from the calibration procedure (Sec¬ 
tion 13.111 . As mentioned in Section [T1 the ionosphere 
contributes most significantly to propagation effects rel¬ 
evant to sky-averaged spectra via absorption, emission, 
and refraction. 

3.5.1. Contribution from ionospheric absorption and 
emission 

The ionospheric absorption, at a given time t, can be 
characterized by the optical depth t{v, 9, (p) (resulting in 
signal loss in the ionosphere Liono{v,9,(j)) = 
and thermal emission of electrons at temperature Tg 
(assumed here to be direction in dependent). There¬ 
fore, based on radiative transfer (|Rvbicki fc LightmanI 
1198511 . the ionosphere-modified sky brightness tempera¬ 
ture, TP{9,(j)), can be expressed as: 

TP{9, P) = Tff9, </,)e-"("’®’« + Tc{l - 

~ Tu{9, (/)){! - t{v, 9, (p)) + Tct{v, 9, (p) 

where the approximation assumes r -C 1. 

We further assume that t{d,9,P), is symmetric in az¬ 
imuth, which due to diurnal gradients in solar irradiation 
is only an approximation. We tested this assumption us¬ 
ing the slant TEC (STEC) data obtained from the Space 
Weather Services (SWS) webpage0, which was computed 
for the MRO location based on GPS satellite data from 
Australian GPS stations. According to these data, the 
symmetry assumption is valid within ~10-20% in the 
North-South direction (slightly better in the East-West 
direction). Moreover, the occurrence of tubular plasma 
over-density structures aligned along Earth’s magnetic 
he ld in the upper i onosphere has recently been reported 
bv iLoi et al.l (1201511 . which occasionally introduce plasma 
density irregularities of the order of 10%. Nevertheless, 
such unusual effects were not taken into account in our 
estimations. 

Therefore, within our approximation, t{i), 9, p) can be 
factorized into a product of a frequency dependent term 
at zenith T2{v) and a frequency and azimuth independent 
geometric term rg{9), which depends only on zenith an¬ 
gle: 


( 2 ) 


http: //www. ips. gov. au/World_Data_Ceiitre/l/1 
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Table 1 

Sunrise and sunset times at the MRO calculated for start, middle and end of every analyzed month. 


Date 

Sunrise 

(AWST) 

Sunrise 

(LST) 

Sunset 

(AWST) 

Sunset 

(LST) 

LST - AWST(hours) 
at local midday 

2014-10-24 

05.48 

07.42 

18.38 

20.30 

1.9 

2014-11-01 

05.37 

07.83 

18.48 

20.91 

2.4 

2014-11-15 

05.23 

08.61 

18.63 

21.99 

3.3 

2014-11-30 

05.17 

09.60 

18.85 

23.25 

4.3 

2014-12-15 

05.22 

10.56 

19.02 

00.33 

5.3 

2014-12-31 

05.35 

11.75 

19.13 

01.51 

6.4 

2015-01-15 

05.53 

12.92 

19.17 

02.53 

7.3 

2015-01-31 

05.75 

14.18 

19.10 

03.51 

8.4 



Frequency [MHz] 

Figure 8. Example of differences calculated between spectra averaged in 0 — 1 h LST ranges and a median spectrum of all 1 h LST spectra 
for the data collected between 2014-10-24 and 2014-11-30. The data are presented in the original 117.2 kHz resolution and some channels 
affected by RFI have been removed for clarity. 



LST [hours] 

Figure 9. Standard deviation of a Gaussian function fitted to dis¬ 
tributions of differences of 1 h LST spectra and the median spec¬ 
trum in several frequency bins. The standard deviation is shown as 
a function of LST. It is stable for 1 h LST bins collected during the 
nighttime and is significantly larger for the 1 h LST bins collected 
at daytime. The seemingly abrupt transition between 97.5 — 102.5 
and 112.5 — 117.5 at LST <8h results from a steep ^ectrum of 
the ionospheric effects at frequencies <110 MHz (Fig. 


/ 100 

t{v, 6, (j)) = T:,{v)rg{9) Ri tioomhz (- rg{6), (5) 

^ I'MHz ' 

where we introduce tiqqmhz = 0.23 x ^ 

0.0023 representing optical depth at zenith at 100 MHz 
and rg{9) = As/Ahjj is a ratio of ray path length As in 
the D-layer to the width Ahjj of the D-layer (Fig. ITO)) . 
The exact equation for rg{9), which was derived from ge¬ 
ometrical considerations of the triangles COU and COB 
(Fig. [ini) using the law of cosines, can be approximated 
by the following formula: 




1 -|- Hd/R f, 

v/cos2(0)+2iJ^^’ 


( 6 ) 


where ij£)=75km is the height of the D-layer above 
the surface of the Earth of radius i?e=6371km, and 
Aho = 30 km is a mean width of the D-layer. Because 
the length of the radio-ray in the ionosphere depends on 
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Figure 10. Schematic of the ionospheric refraction in the F-layer 
and absorption/emission in the D layer. Due to the ionospheric 
refraction the radio-ray intercepted by the antenna at zenith angle 
6 originates from a slightly lower part of the sky at zenith angle 

e + 56{u,e). 


its path, the value of rg {9) factor is also slightly modified 
by ionospheric refraction, but we estimated this effect to 
be negligible as the path length change is <0.1% even for 
relatively high refraction effects ~larcmin at 100 MHz 
and 45° zenith angle. 

The antenna pattern of a conical log spiral antenna 
should, by design, be frequency independent. Agree¬ 
ment between expected antenna properties and the real 
antenna was confirmed by comparison of calibrated sig¬ 
nal and predictions of the sky model integrated with the 
FEKO-simulated antenna pattern and by the measured 
reflection coefficient of the conical log spiral antenna in 
the 70 —200 MHz frequency range (Fig.[^. Therefore, in 
order to simplify our derivation, we further assume a fre¬ 
quency independent antenna pattern, Pv{0, 4>) = P{9, 4>). 

After substituting equation [5] into equation |4] and 
then equation 2] into equation |31 we can express the 
ionosphere-affected antenna temperature Ta„t(iz) as a 
sum of “no-ionosphere” antenna temperature 
and a ionospheric perturbation term: 


where substituted according to equation [2] 

and rg/ is the P^{9, (j))T,y{9, (/))-weighted average of rg{9). 
After defining rg = ^ P^{9,(j))rg{9)dV,, equation 0 

can be expressed as: 


/ 100 \2 

Tant{v) ~ IQOMHzTe (-) rg 

/ 100 \2 
-TlOQMHzi -) 

Based on the FEKO simulated antenna pattern, the 
value of the frequency independent geometric factor rg 
was found to be r™”® Ril.4 and, as expected, almost 
constant at frequencies above 75 MHz. For comparison 
we have also calculated Ri 1.7 for a Gaussian sym¬ 

metrical beam ea;p(—(0/52°)^)), which we found to 
be a good and smooth approximation of the FEKO sim¬ 
ulated pattern of the conical log spiral antenna (Fig [3]). 
Similarly, the value of frequency independent geometric 
factor rg/ for the cone was found to be rg/ r; 1.48. 

The contribution from ionospheric absorption and 
emission to the difference between the spectrum T^{v) 
collected on the day and a reference spectrum T” {v ), 
both collected in the same LST range but at different 
ionospheric conditions characterized by two pairs of pa¬ 
rameters and (t[oomhz>> can be ex¬ 

pressed by the following formula: 


II 

<1 

T^ntiP 

)-TLt{’/) = 


(■^lOOMHz ~ ^lOOMHz ) ( 

100 N 

^TZt{//)rg/+ 


//MHz/ 

(10) 

{^n rpn r rpr\ 

VlOOMHz-^e ^lOOMHz-^e } 

i( )V.. 

vmhz 7 



If we approximate the sky temperature with a power- 
law T™Ai') ~ Tinn(10 0/i^MHz)°. where a ^ 2.6 

(iRoerers fc Bowm;ml 12008( 1 and Tioo is the sky temper¬ 
ature at 100 MHz, we obtain the expected frequency de¬ 
pendence of the difference between two spectra collected 
at the same LST expressed as a sum of two power-laws 
~ \jv^ and ~ l/iz^+“: 


A“(^) 


rg^Tr 


/ 100 ' 
^ r'MHz - 


- Tg/ArTioO ^ 


100 \2-Ha 

? 

( 11 ) 


I'M Hz/ 


/ 100 \2 1 

Tanti//) +P00MHz( -) X 

^ I/MHz ' t'Q 

( [ P,{9,^)Terg{9)dn- [ P,{e,cj))T,{e,(l))rg{9)dn). 
\ J 4:77 j 477 / 

(7) 

The second integral can be expressed based on the 
mean value theorem for integrals as: 


p,{o,4,)T,{e,/i>)rg{e)dn = 

G) J47I- 

V ^ 0)T.(0, = rg/ T“,(zz), 


( 8 ) 


0 J 477 


where we define = (r” -H"oomhz) and Att = 

(-n rpn _ r rjirS 

VlQOMHz^e ‘lOOMHz/e )■ 

Because the values of rg and rg/ are very similar we fur¬ 
ther approximated rg/ k. rg k, 1.4, which simplifies the 
following analysis without impacting the frequency de¬ 
pendence. Nevertheless, this approximation has a small 
(<5%) impact on the measured electron temperatures. 

3.5.2. Contribution from ionospheric refraction 

Ionospheric refraction is relevant to global EoR exper¬ 
iments because it distorts the sky observed by a ground- 
based instrument such that the ionosphere acts as a huge 
slightly dc-magnifying lens. Therefore, from the antenna 
perspective, a certain pointing {9,(j)) within the antenna 
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beam gets deflected by the ionosphere such that the an¬ 
tenna intercepts the radiation originating from a slightly 
lower sky position {6 + S 9(i/,d ), 4>) (Fig. flOl) . 

As described in Section ITTTl the typica l mean offset ex- 
pected from the theoretical formula by iBailevI (jl948fl is 
S9q ~ O.larcmin at 100 MHz and zenith angle 45° and 
varies proportionally to variations in electron density. 
Similarly, based on the data fro m the MWA telesc ope 
(co-located with our instrument) lArora et al.l (|2015D re¬ 
port nighttime offsets of the source positions of the order 
of tenths of arcmin at 100 MHz and scaling as 1/z/^. 

Therefore, based on the TEC values computed from 
the GPS data for the MRO (the same data from the SWS 
webpage as in Sec. 13.5.IL we estimate that the intra¬ 
day variations of the reference offset S9o within a month 
should typically be ^20-30% and can be higher during 
the periods of higher solar activity. Bailey’s formula for 
the deflection angle can be expressed as a product of 
frequency and zenith angle dependent terms as: 


69{iy,9) = 59{v)89{9) ^ 59o{im/vMHzf (12) 

where 59q ~ O.larcmin. The exact form of fre- 
qu ency independent geometric fact or f{9) can be found 
in (IBailevI[l948l : iDatta et al.ll2015ll . After including re¬ 
fraction, equation [2] can be re-written as: 


TZti^) = ^ [ P.(d»r,(0 + <5%,0),0)dH, (13) 

where S9{i/,9) can be calculated from equation [121 
In order to estimate the magnitude of refraction, we 
implemented the above equation in the modeling soft¬ 
ware and calculated sky-averaged spectra for the “non- 
ionospheric” case for several refraction 

cases represented by different values of S9o 

ranging from O.larcmin (typical nighttime) to 5arcmin 
(extreme daytime case). Next, we calculated differ¬ 
ences between these sky-averaged spectra. The dif¬ 
ferences between spectra affected by typical nighttime 
refraction {S9o = O.larcmin) and the hypothetical “non- 
ionospheric” scenario, compared with similar differences 
for absorption/emission effects are shown in Figure [TTJ 
The contribution of nighttime refraction to the sky- 
averaged spectrum above 80 MHz is ^1K, which is much 
less than the expected contribution of absorption and 
emission. 

We also calculated differences between sky-averaged 
spectra affected by different refraction conditions {S9o = 
0.1 — 0.3 arcmin) and found out that even relatively large 
(~0.1 arcmin) variations of S9o around a mean value lead 
to a very small (<1 K) changes of the sky-averaged spec¬ 
tra. However, the same tests for daytime conditions 
(with S9o in the range 1 — 5 arcmin) show that during 
daytime contribution of the ionospheric refraction may 
be of the order of tens of Kelvins. 

In order to derive the frequency dependence of the 
ionospheric refraction, equation 1131 can be approximated 


by the first two terms of its Taylor expansion: 

= PA0, mio. 8^)dn+ 

S9ii^)^ f PA0,^)^^^d^^69i9)dn, 

M) CLU 

where the first term is the antenna temperature T'“((i/) 
defined in equation [2] as measured without (or outside) 
the ionosphere and we expressed S9{h', 9) as a product 
of frequency and angle dependent terms (Eq. [12]). It 
is worth noting that the second term vanishes when 
dT^{9, (p)/d9 = 0, meaning that refraction effects van¬ 
ish in a hypothetical “constant sky” scenario (i.e. when 
T^{9,<j)) is constant). 

Assuming, as before, a power-law dependence of 
the sky temperature T^{9,(p) = T'ioo(d, ^!')(100/i/Mffz)“ 
equation [T71 can be approximated by: 

TZP{v)^TZt{i^)+59{v)( -) hr. (15) 

VmHz ' 


where we define a frequency independent integral IdT as 

IdT = ^ [ PZ0,<l^)^^^^^j^69{9)drt. (16) 

Po Jatt 

Hence, the refraction induced difference be¬ 
tween the spectrum collected on the day 

and a reference spectrum T”®l’'’'’(i/), both collected at 
the same LST but at different ionospheric conditions re¬ 
sulting in different refraction angles dd" and S9^, can be 
expressed by the following formula: 


. /I no \ 

AZP{>y) = {69Z’^)-S9Z^))( -) IdT 

^ VMHz ' 

/ 100 \2+a 

^(S9^-S9^o)( -) IdT, 

^ VmHz ' 


(17) 


where we used equation[T2|to expand frequency depen¬ 
dence of 59'^{v) and S9'~{i'). 

The exponent (2 -|- a) in equation [T7| is exactly the 
same as the exponent in one of the components in equa- 
tion[TT| Hence, we conclude that the difference Anr{v) = 
of two spcctra n and r collected at the 
same LST, due to ionospheric absorption, emission and 
refraction, can be described by a sum of power-laws: 


Anr{v) = A{1QQ/vmHz)^ + BZQQ/vmHz)^ Tant{v) 

^ A{imlvMHzf + H'(100/zzm44.)2+“. 

(18) 


3.6. Comparison of expected ionospheric contributions 

In order to compare contributions of the ionospheric 
effects, they were incorporated into the signal modeling 
software described in S15. Ionospheric absorption and 
emission were included according to equations [3] [4| and [5] 
whilst refraction was included according to equation 1131 
Because the antenna gain values were calculated by 
the FEKO software in discrete directions (every 1°), 
for the sake of this comparison we used the Gaussian 
beam described above. The sky model was integrated 
with the Gaussian beam in the absence of any iono¬ 
spheric effects {TZti^)), in the presence of absorption 
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Figure 11. Example simulations of differences of sky-averaged 
spectra affected and unaffected by ionospheric effects generated for 
24 hours of local time (different curves represent different LSTs). 
The dashed red lines rep res ent the contribution from ionospheric 
refraction (equations 1131 and 1 12 II and solid blue lines represent the 
overall change of received power due to ionospheric absorption and 
emission (equations [S] [HIS] with = 0.01 dB). 


and emission and in the presence of refraction 

only The expected contribution from absorp¬ 

tion/emission was calculated as ^modeii^) = ^ant(^) ~ 
T™t i'^) the expected contribution from refraction as 

(Fig. mi. 

Similarly, we also calculated differences between typ¬ 
ical and perturbed parameters (by amounts expected 
from the TEC variations at the MRO) and the results 
also support the assertion that absorption and emission 
will dominate at these fre quencies, which is also in agree¬ 
ment with simulations bv iVedantham et al.l (|2014f) . Dur¬ 
ing the nighttime, the contribution from refraction at 
80 MHz is almost negligible (below 1K) in comparison to 
absorption and emission when mean offset and its vari¬ 
ations at 100 MHz and zenith angle 45° are ~0.1 arcmin 
and slightly larger (<5K) for ^0.5-1 arcmin. Thus, we 
expect that at least during the nighttime, fitting of dif¬ 
ferences of the spectra with equation [TH] will result in 
values of parameters A and B being dominated by ab¬ 
sorption and emission with a negligible contribution from 
refraction. 


3.7. Fitting ionospheric parameters 

Assuming that the differences between each 1 h spec¬ 
trum and the median spectrum (Fig. |8]) are due to vari¬ 
ations of the sky signal caused by the ionospheric effects 
described in Section 1331 we fitted the difference spectra 
with equation [181 Example fitting results for the LST 
0 — Ih bin collected in Oct and Nov 2014 are shown in 
Figure [12] and Table [2] 

3.7.1. Derivation of electron temperature 

We plotted the fitted A values as a function of B values 
(eq. [TH]) ; an example of their dependence for the 0 — 1 h 
LST bin is shown in Figure [131 which suggests a linear 
correlation between these parameters (at least at some 
LSTs). The fitting procedure was repeated for all LST 
ranges and the results also exhibit a very similar linear 
correlation, particularly for the nighttime data. Based 
on equation HU the fitted parameters A and B can be 
interpreted as a combination of frequency independent 


Table 2 

Parameters A and B, thei r er rors, and reduced resu lting from 
fitting equation 1181 to the data in Figure [TJ] 


Date 

A 

AA 

B 

AB 

X^/NDF 

20141129 

0.0114 

0.0007 

3.59 

0.41 

60.7 / 24 

20141127 

0.0293 

0.0009 

16.21 

0.48 

207.7 / 24 

20141123 

0.0083 

0.0007 

5.38 

0.39 

32.2 / 24 

20141121 

-0.0017 

0.0007 

-2.89 

0.41 

17.1 / 24 

20141120 

0.0091 

0.0007 

5.73 

0.40 

32.0 / 24 

20141118 

0.0001 

0.0006 

1.40 

0.37 

14.0 / 24 

20141116 

-0.0100 

0.0008 

-4.51 

0.46 

151.6 / 24 

20141114 

-0.0212 

0.0009 

-5.80 

0.47 

186.9 / 24 

20141112 

-0.0131 

0.0007 

-5.56 

0.43 

79.5 / 24 

20141110 

0.0026 

0.0006 

6.51 

0.39 

27.0 / 24 

20141108 

-0.0034 

0.0006 

-3.47 

0.42 

20.7 / 24 

20141106 

0.0033 

0.0006 

2.65 

0.40 

34.3 / 24 

20141103 

-0.0006 

0.0006 

-0.99 

0.40 

15.5 / 24 

20141101 

-0.0047 

0.0006 

-3.81 

0.39 

39.1 / 24 

20141031 

-0.0193 

0.0008 

-10.63 

0.45 

158.6 / 24 

20141029 

0.0045 

0.0008 

4.54 

0.54 

21.1 / 22 

20141026 

-0.0102 

0.0008 

-5.35 

0.41 

67.1 / 24 

20141024 

0.0006 

0.0006 

-0.29 

0.36 

18.9 / 24 


parameters 

B = i'^lOOMHz ~ '^100MHz)'''gh 

^ = i'^lOOMHz'^e ~ ''"lOOMHz^e )’'g 
~ F'e(r"gQJ^^^J, — Ri TeB, 

where we assumed that for a given 1 h LST range the 
electron temperature is approxi mately constant T" r; 
Tf R! Te and rgt r; rg (see Sec. I3.5.1|l . Hence, in this 
approximation A k, T^B and therefore the slope fitted 
to the data in Figure [13] can be interpreted as electron 
temperature for a given LST range. 

As further results will show, the constant Tg assump¬ 
tion is valid for the exclusively nighttime 1 h LST bins, 
but it is invalid for other LST bins. Because LST drifts 
with respect to local time by r2 h/month, a 1 h LST bin 
effectively covers r3 h in local time (the width of the 1 h 
LST bin plus 2 h of the drift between LST and local time) 
over the period of a month (3.5 h in the case of Oct/Nov 
dataset). The assumption is invalid particularly for the 
LST bins near sunrise/sunset when electron temperature 
changes rapidly (but smoothly) or for the daytime bins 
when electron temperatures can vary between days. The 
results for these LST bins therefore represent an aver¬ 
age electron temperature in a 3h (3.5 h in the case of 
Oct/Nov data) local time bin (as indicated later by the 
bars on the time axis). 

The good correlation between the fitted parameters A 
and B indicates that ionospheric absorption and emission 
dominate variations of the sky-averaged spectra whilst 
ionospheric refraction is less significant, which agrees 
with our initial estimates in Section 13.61 As mentioned 
above, a very good linear correlation was typically ob¬ 
tained for nighttime data, but daytime data show much 
larger scatter (Fig. HU d,e,f,g). There may be multiple 
reasons for this, such as non-perfect excision of data af¬ 
fected by the Sun or that other effects (i.e. ionospheric 
refraction) play a more significant role during the day¬ 
time. Nevertheless, we performed the same analysis for 
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80 100 120 140 160 180 200 


Frequency [MHz] 

Figure 12. Example results of fitting equation ^8]to differences of 
0 — Ih LST spectra (as in Figure [8| after averaging of 40 frequency 
channels) and a reference median spectrum collected in Oct and 
Nov 2014. The spectra are offset by 50 K for clarity, but they 
all lie around zero. Large errors are usually due to residual RFI 
(for example due to the FM-band in the case of two points in 
90 — 100 MHz band on the curve at f^lOOK). The values of fitted 
parameters A and B, their errors, and reduced ^re listed in 
Table [2] 


all 1 h LST bins for three calendar months and derived 
Te as a function of LST (Fig.[T^. The dependence of Tg 
on local time will be presented in Section |4l 

3 . 7 . 2 . Derivation of variations in ionospheric absorption 

According to equation [19] the parameter B, fitted to 
differences of 1 h LST spectra, can be interpreted as 
B = {twomhz - TiqoMHzyg- Thus, its distribution 
represents variations in the ionospheric optical depth at 
100MHz (tioqmhz)- For every Ih LST bin we created 
a histogram of fitted parameter B. An example distri¬ 
bution of B for LST bin 0 — Ih in Oct/Nov 2014 data 
(distribution of X-axis in Fig.fT^l is shown in FigurefTKl 
Because Figure [T5l is a distribution of differences of opti¬ 
cal depths at 100 MHz h z ~'^womh z)^ Ts standard 

deviation is larger by than the standard deviation of 
the optical depth itself ar = o'(riooMifz)- Therefore, in 
order to calculate the standard deviation of tiqqmhz at 
zenith, the calculated standard deviation of parameter B 
was divided by a factor V2rg « 2, where rg Ri 1.4 is the 
value of the geometric factor calculated for the conical 
log spiral antenna (Sec. I3.5.T]) . 

The standard deviation of optical depth tiqqmhz as a 
function of local time calculated for all 1 h LST bins and 
all four months is shown in Figure 1161 The observed 
nighttime optical depth variations are typically CTt- ^ 
0.005, which corresponds to ^0.002 at 150 MHz . very 
similar to values measured by [Rogers et al.l (j2015ll . Al¬ 
though our value, CTt- ~ 0.005, is slightly larger than typ¬ 
ically cited nighttime optical depth at zenith tiqqmhz 
0.0023 (~0.01dB), it is well within a range of values 
(~ 0.0023 — 0.023) cited in the literature (Sec. I3.5.1|) . 
Moreover, optical depth and its variations may also be 
season and solar cycle dependent. 

3 . 7 . 3 . Zenith absorption at 100 MHz 

In order to measure the optical depth at zenith 
at 100MHz {tiqqmhz) we used a riometer-like ap¬ 
proach. Based on multi-day observations a quiet 
day curve (QDC) of the maximum observed sky tem¬ 
perature (at a state of very small residual absorp¬ 
tion) is determined and rela tive absorption is measured 
with respect to the QDC (iKrishnaswamv et al.l 119851 : 
iBelikovich fc Benediktovlll969fl . 

Due to low-level RFI in the FM band (88—108 MHz) we 
determined the optical depth at the RFI free frequency 
of 81 MHz and scaled the result to 100 MHz by a fac¬ 
tor of (81/100)^. The calibrated antenna temperature 
at 81MHz (r8,lf«-) was calculated as a mean value of 
Tant in the 80 — 82 MHz sub-band. In order to calcu¬ 
late absorption in the selected band, we first determined 
a complete (for nighttime and daytime) QDC based on 
1 h LST spectra from all the data collected between Oct 
2014 and Jan 2015. The same quality selection crite¬ 
ria as described in Section 13.71 were applied in order to 
excise data affected by any kind of variability (mainly 
solar or lightning). In order to determine the QDC value 
for a given LST bin we used a sample of antenna tem¬ 
peratures from quality-accepted 1 h LST spec¬ 

tra. If for a particular LST the sample contained spectra 
collected during nighttime and daytime, only the night¬ 
time bins (solar elevation < —2°) were used. The QDC 
value TZTiLST) = calculated for all 
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(a) LST = 0-1 hours 








(g) LST = 20-21 hours 



Figure 13. An example of correlation between fitted para meters A and B for several 1 h LST bins for the data collected in Oct/Nov 2014. 
Parameters A and B can be interpreted as in equation ca and fitted values of Te are shown in each figure. An approximate range in local 
time can be obtained by subtracting P=J3.5hours. The small number of data points in figures c) and g) is a result of exclusion of Ih LST 
bins including sunrise and sunset or containing nighttime and daytime spectra. 
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Figure 14. Electron temperature as a function of LST for 
Oct/Nov 2014, Dec 2014, and Jan 2015 (see Table^for translation 
between LST and local time). 



Figure 17. Quiet day curve (QDC) derived as maximum value of 
antenna temperature where is a mean 

value of Tant in 80 - 82 MHz Mn. 




Figure 15. Example distribution of fitted values of parameter 
^ = (’'lOOM/fj - the LST bin 0 - 1 hours from 

Oct/Nov 2014 data (X-axis values from Fig. I13k b 



Figure 16. Standard deviation of zenith optical depth at 100 MHz 
as a function of local time derived for Oct/Nov 2014, Dec 2014, and 
Jan 2015 data. 


Ih LST bins (including the ones collected only during 
daytime) is shown in Figure [TTl We further assume that 
the QDC curve, T^^{LST), represents antenna tem¬ 
perature measured at a time of very small residual iono¬ 
spheric absorption, and we determine relative absorption 
with respect to this QDC curve. Due to limited amount 
of data and the fact that the sta te of minimum absorption 
occur s in the pre-dawn hours (IBelikovich fc Benedikto^ 
I1969I1 . this assumption is valid for only some nighttime 
LST bins and invalid for LST bins collected during the 
daytime only. 

Based on equation 0] we can define the difference be¬ 
tween the maximum measured antenna temperature at 
81MHz TJ^^{LST) (at negligible absorption) and any 
other measurement Tant{LST) as 

Adc = TZT {LST) - (LST) = 

Ti,,,LST){TZT{LST)-T,{LST)), ^ ^ 


Figure 18. Optical depth at 100 MHz as a function of LST derived 
according to variations of antenna temperature (mean 

value of Tant in 80 — 82 MHz bin) with respect to its maximum 
value T^^^{LST). The results for Ih LST bins collected only 
during daytime were excluded from the analysis. The errors are 
~ 0.005 and therefore c onsist ent with the optical depth variations 
(ctt) derived in Section 13.7.21 using a slightly different method. 


and therefore, the optical depth at any given frequency 
and LST t{v,LST) can be calculated from a difference 
between QDC and other spectra as: 


r(iz, LST) 


TZrjLST) - T^^,{LST) 
T^^-{LST)-T,{LST) ’ 


( 21 ) 


where Te{LST) is the electron temperature at a given 
LST and was calculated according to the parameteriza¬ 
tion of the data in Figure[T4l Equation [2T] indicates that 
the accuracy of r(zz, LST) strongly depends on the accu¬ 
racy of Tf,[LST), which is very good for the nighttime 
data, but varies significantly for daytime LSTs. There¬ 
fore, we only determined LST) for LSTs which could 
be collected during nighttime, which excluded LST bins 
14-21 h. 

In order to estimate a commonly cited optical depth 
TiooMHziLST) at zenith we divided t{81MHz, LST) by 
the geometric factor r™”® r; 1.4 and subsequently multi¬ 
plied by a factor (81/100)^ « 0.66 (eq. [5]). For each LST 
bin, values of tiqqmhz were calculated from de¬ 

rived from all 1 h spectra and their mean and standard 
deviation as a function of LST are shown in Figure fTBl 

In order to cross-check our result, we derived tiqqmhz 
from other frequency sub-bands (92 — 94 and 102 — 
104 MHz) and the results were very similar and within 
the measurement error. The errors in Figure [THl are 
R! 0.005 and therefore consistent with the standard devi¬ 
ations of optical depth (cr,-) derived in Section [3.7.2l using 
a slightly different method. 
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The measured values of nighttime optical depth 
tioomhz ~ 0.01 (loss of 0.05 dB) at zenith are about 
a factor of 4 larger than the typically cited value 
-0.0023 (-0.01 dB), but within a 0.0023 - 0.023 (0.01 - 
0.1 dB) range of values c ited in literature for nighttime 
(lEvans fc Hagfor.^ Il968t iDaviesI Il99(11 iThomoson et al.1 
l2007fl . Moreover, the measurements performed with tra¬ 
ditional riometers did not take emission into account 
(Te{LST) = 0 in equation [2T|) . Therefore, our values 
of optical depth can be slightly higher than those cited 
in the literature, but we find the presented method more 
accurate, especially for measuring very small nighttime 
optical depths. For comparison we also calculated opti¬ 
cal depths according to equation [21] with Te{LST) = 0 
(r(i.,L5T) = (T^^riLST) - T^^,{LST))/TZr{LST)) 
and it shifted the optical depths from —0.01 to —0.007. 
Since nighttime absorption can also be season and so¬ 
lar cycle dependent, we consider our measurements to 
be in satisfying agreement with other measurements at 
the same loc ation, but different time of year and smaller 
data sample ([Rogers et al.ll20lll . Moreover, an analysis 
based on a larger data sample should allow us to further 
verify and refine our results. 

The presented absorption analysis can be considered as 
a proof of principle because it was based on a relatively 
small data sample (3 months worth of data), which is not 
sufficient for precise determination of a quiet day curve 
(QDC) in all LST bins. Present analysis was strictly 
valid for a few nighttime LST bins. The complete QDC, 
built from the one year of data, will enable us to cor¬ 
rectly determine relative absorption in all LST bins. A 
complete QDC may be shifted towards higher values and 
result in slightly higher values of optical depth, but this 
should mostly affect LST bins presently observed during 
the daytime. We consider our present nighttime mea¬ 
surements of optical depth reliable and attribute their 
slightly higher than expected values primarily to time of 
the year and solar activity dependence. 


4. COMPARISON WITH IONOSPHERE MODEL 
PREDICTIONS 

We compared our measured electron temperatures 
with predictions of the Inter national Reference Iono¬ 
sphere 2012 (IRI-2012) model ([Bilitza et al.ll20i^ . We 
initially derived the dependence of Tg on LST (Fig. [T4|). 
Therefore, in order to derive the dependence of Tg on lo¬ 
cal time for each calendar month, we converted LST to 
local time by applying the appropriate shift calculated for 
the middle of the given calendar month (Tab. |T|). The 
measured electron temperatures as a function of local 
time for Oct/Nov 2014, Dec 2014, and Jan 2015, over¬ 
plotted with the predictions of the IRI-2012 model for 
2014-11-15 weighted by optical depth and at four dif¬ 
ferent heights (60, 140, 500 and 1000km), are shown in 
Figure [T9| The overall trend of the measured and model 
curves is similar as they both show a rise in electron tem¬ 
perature at sunrise and fall at sunset, but the values of 
measured and predicted temperatures differ (especially 
during the daytime). 

The weighted electron temperature (Tg) was derived 
from the results of the IRI-2012 model according to the 


following equation: 

^ rg(h)jVg(fe)i.g(h)dfe ^ 

jhmax ]Vg(/l)z/g(/l)(i/l 

where I'dh) is a parameteriza t ion of collision rate data 
(Figure 7 by lAggarwal et al.l l|1979| l) as a function of 
height h, h^ax = 2000 km is the maximum height of the 
data generated with the IRI-2012 model, and Ne{h) is 
electron density as a function of height obtained from the 
IRI-2012 model. We found, however, that for h < 80 km 
the electron densities according to the IRI-2012 model 
were underestimated (typically Ng = 0e~lm^ for night¬ 
time). Therefore, nighttime Tg was originally biased to¬ 
wards higher values by contribution from higher iono¬ 
spheric layers, because there was no (or insignificant) 
contribution from the D-layer. Hence, whenever Ne(h) 
was smaller than typical nighttime resi dual electron den¬ 
sity in the D-layer. A7i = 10^e~/m^ (lEvans fc HagforsI 
1 1 968t iDaviesll 1 99(i iThompson et al.ll2007ll . we artificially 
set Nf. to Nq. This correction enabled a very good agree¬ 
ment between our nighttime measurements of electron 
temperature and weighted electron temperature derived 
from the IRI-2012 model. 

The measured nighttime electron temperature is ap¬ 
proximately constant (~ 470 K) and lies between model 
curves generated for heights 60 — 140 km (D and E- 
layers). This is the region expected to primarily con¬ 
tribute to the nighttime ionospheric absorption, whilst 
higher layers are not expected to contribute significantly 
to absorption during nighttime, and therefore to our 
measured value of the electron temperature. We consider 
our nighttime measurements as the most reliable because 
they are not affected by solar activity and also our results 
show that the assumption about the dominant role of ab¬ 
sorption holds best for the nighttime data (Sec. 13.7.11) . 
The nighttime data are also particularly interesting from 
the perspective of the global EoR experiments. 

The electron temperatures measured during the day¬ 
time, reaching as high as ~3000K, are significantly 
larger than the weighted electron temperatures derived 
from the IRI-2012 model, which are Tg ^1000 K. The 
disagreement may be a combination of several factors. 
Firstly, our method may not be applicable to daytime 
data because of solar contributions affecting our data 
analysis (solar activity that was not excised, and move¬ 
ment of the Sun within the beam pattern). Secondly, 
the larger contribution of F-layer refraction is not taken 
into account in the presented analysis. Thirdly, daytime 
measurements lie mostly between model predictions for 
heights 500—1000 km (Fig.[T9|), which might be a result of 
significant daytime absorption in the F-layer. Thus, lack 
of exact knowledge of collision rates at heights h>500 km 
(approximated in our calculation of Tg by 100 Hz ) may 
lead to an underestimation of Tg. 

Therefore, we conclude that our measurements agree 
with the IRI-2012 model predictio ns (particularly at 
nighttime) and othe r measurements (jZhang et al.ll200'^ 
iPawsev et al.l Il951t) within the applicability of our 
method. The agreement can be tested better in the 
future when changing quiet-Sun conditions, un-excised 
solar activity remaining in the data sample, and iono¬ 
spheric refraction are better accounted for in the data 
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Figure 19. Comparison of measured electron temperature (data 
points with error bars) w ith predictions of In ternational Reference 
Ionosphere 2012 model tlBilitza et al.l [20141^ weighted by optical 
depth (black solid line) and at four heights: 60, 140, 500 and 
1000 km (different dashed lines). The bars on the time axis corre¬ 
spond to ^3 h local time coverage of a 1 h LST bin over a period of 
1 month (3h results from Ih LST bin itself plus 5^2h/month drift 
of LST with respect to local time). 



Figure 20. Calibrated sky temperature at 80 MHz as 

a function of LST overplotted with its parameterization by a 
Fourier series (eq. I23II and the maximum («"") and mini¬ 
mum values of the fitted Fourier series. We note that 

1.9 — 3.9 h LST appears to be the best range in which to integrate 
data down for the global EoR experiment. 


analysis. 


5. STABILITY TEST OF SKY DATA 

According to the dependency of with local 

time (Fig. [201), we identified that was minimal 

at around LST = 2.9 hours and not at Ri5.8h when the 
Galactic Center was at its lowest point below the horizon. 
Therefore, we decided that the optimal conditions to se¬ 
lect spectra for averaging are at LST^pt ± ALSTopt = 
2.9 ± 1 hour and only during the nighttime, which re¬ 
sults in <2h of data per day to be used for the global 
EoR analysis (as some data might be excised by RFI or 
variability criteria). 

In order to test whether the sky spectra will integrate 
down to the required precision, we calculated the stan¬ 
dard deviation of Tant from samples of increasing num¬ 
bers of integrations. The standard deviation (cts) was cal¬ 
culated from a sample of spectra in the reduced and cal¬ 
ibrated dynamic spectrum collected in the period 2014- 
10-24 to 2014-12-31 (in 35sec resolution). We could not 
include the data after 2015-01-01 because the optimal 
LST range (1.9 - 3.9 h) began to enter daytime. Then, for 
several selected frequency channels, as was plotted as a 
function of integration time Tint, related to the number of 
sample integrations Ng as Tint = Ns x 35 sec. Such a test 
allowed us to verify that after a sufficiently long integra¬ 
tion time, CTs converges and remains approximately con¬ 
stant with an increasing number of sample integrations. 
Therefore, averaging an increasing number of spectra re¬ 
duces the statistical error as expected from equation [T] 
by a factor ^ l/^/N^ (Fig. [^, which is very encour¬ 
aging and convinces us that ionospheric or instrumental 
effects should not prevent us from integrating down to 
the required precision. 

5.1. Expectations based on system temperature 

In order to verify our understanding of the as depen¬ 
dence on the integration time we generated a model dy¬ 
namic spectrum with separate time and frequency com¬ 
ponents and use this model to generate simulated obser¬ 
vations. Firstly, we parameterized the sky temperature 
dependence on LST [Tant{LST)) with a Fourier series: 

6 6 

F{LST) = flQ + aiSin{2'Kip) + biCos(2Trip), (23) 

i=l i=l 



Figure 21. The standard deviation of Tant (cts) divided by the 
square root of the number of samples {Ns = Tint/3b sec) ^ a func¬ 
tion of integration time Tint for several selected frequency chan¬ 
nels. Because cts converges and remains approximately constant 
when Ns is sufficiently large, the standard error of the mean de¬ 
creases as expected from equation ^ by ^ factor ~ 1/ y/Ns (blue 
line). Therefore, the required precision may be reached after suffi¬ 
ciently long integration time. The test was based on 2014-10-24 - 
2014-12-31 data. The data collected after 2015-01 could not be in¬ 
cluded because the optimal LST range (1.9 - 3.9 h) began to enter 
daytime. 



Figure 22. Sky spectra at maximum (T^°J(i/)) at LSTf^l7.8h 
and minimum {T^nt^{i^)) at LST~1.9h parameterized with sums 
of two power laws. The sum of two power laws yielded a slightly 
better fit than a single power law with fit residuals within P=i20K 
for “hot” and p^IOK for “cold” sky over most of the bandwidth. 


where p is (LST/24). The data with a fitted Fourier 
series are shown in Figure [201 where we also defined the 
maximum and the minimum values of 

the fitted Fourier series. We define a normalized Fourier 
series as: 


F{LST) 


F{LST) - 


7^80MH2 _ j^SOMHz ■ 
max min 


(24) 
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Next we fitted a sum of two power laws in frequency 
to the maximum LSTRil7.8h) and mini¬ 
mum at LST~1.9 h) sky temperature (Fig. [22]). 

We used two power laws because the model yielded a 
slightly better fit than a single power law, with fit resid¬ 
uals within Ri20K for “hot” and RilOK for “cold” sky 
over most of the bandwidth. The antenna temperature 
was generated for the same LSTs as the sky data 
by adding a difference of sky spectra ~ 

modulated by F[LST) to the “cold” sky spectrum ac¬ 
cording to the formula: 

(25) 


This model represents the “real” sky within ±20 K, 
which is sufficient to estimate the errors expected from 
cosmic noise as a function of increasing integration time. 
The statistical error on the simulated observations was 
calculated according the same formula which we use to 
calculate the error on calibrated antenna temperature 
derived from the reference source (equation 13 in S15, 
where the ratio of antenna to reference data was ap¬ 
proximated by one and a very small contribution from 
the error on ambient temperature measurement was ne¬ 
glected): 






(26) 


where B is the frequency bin of Rill7.2kHz, Trcv ~180K 
is the receiver noise, Tant=^^ sec is the integration time 
on the antenna, and r,.e/=12.5sec is the average inte¬ 
gration time on the reference (both from the reduced 
dynamic spectrum). 

The dependence of the standard deviation erf®" on inte¬ 
gration time resulting from the same procedure as previ¬ 
ously applied to our data is shown in Figure |23l The val¬ 
ues of erf®" converge to slightly lower values than corre¬ 
sponding curves in Figure [211 which might be attributed 
to variations introduced by the ionosphere (Fig. 12]) not 
taken into account in this simplified simulation proce¬ 
dure. In order to verify this hypothesis we calculated 
o-tot = ((o-f®")^ ± erf^no)^^^: where aiono was taken from 
Figure [5] at the closest frequency to those shown in Fig¬ 
ures |2T] and |23l Comparisons of as/\/Ws and atotj'/^s 
are shown in Figure |24l Although atot/\/Nl is slightly 
above Ug/^/Wl the overall agreement between the two is 
very good. Assuming variations shown in Figures [5] and [HI 
result mainly from ionospheric variability (at least dur¬ 
ing the nighttime), which we show in SectionsjSjandS] we 
conclude that we understand the magnitude of sky signal 
variations. Furthermore, it turns out that these varia¬ 
tions are dominated (especially during the nighttime) by 
ionospheric absorption and emission. 

The standard deviation calculated from the large data 
sample (~98days) converges to a constant value, which 
proves that, in principle, the additional variability in¬ 
troduced by the ionosphere does not impede integrat¬ 
ing down to the required precision with ground-based 
instruments. Therefore, because ionospheric stochastic 
error can be averaged to a mean ionospheric contribu¬ 
tion, the latter can be identified and removed from the 



Figure 23. Standard error of the mean simulated antenna tem¬ 
perature {(7s/y/Ns) as a function of integration time Tint (number 
of integrations in the sample can be calculated as Ns = Tint/^^ sec) 
for the same set of freq uen cy channels as for the same test applied 
to the real data (Fig. I21II . The s imu lated data are slightly be¬ 
low corresponding curves in Figure [211 which might be attributed 
to variations introduced by the ionosphere (Fig. O not taken into 
account in the simplified simulation procedure. The ripples on 
timescale of ^2 h are due to small changes in sky temperature in 
the optimal LST range (1.9 - 3.9 h). 



Figure 24. Standard error of the mean antenna temperature 
{cTs/\/N s) as a function of integration time Tint for sky data (black 
curve is the same as in Fig. 12111 compared with simulated data 
with additional ionospheric variability (blue curve) calculated as 
ertot/VN^ = where Uiono was taken 

from Figure[9]at the closest frequency to the presented frequencies 
(80, 90 and 125 MHz). 


sky-averaged spectrum (on a daily or hourly basis) in or¬ 
der to minimize the extra variance which it introduces. 
One possibility, for instance, is to identify the ionospheric 
contribution to the sky-averaged signal on an hourly ba¬ 
sis (as in Fig. IE|)- Because daily contributions to the 
total integration time are going to be of the order of a 
few hours (here we tested the LST range 1.9 - 3.9h), 
the ionospheric contribution can be identified by fitting 
equationjTHjto parameterize the differences (as in Fig. US]) 
between I h LST spectra and a median spectrum of this 
I h LST bin. Then, such a parameterization would be 
subtracted from the original 1 h spectrum before averag¬ 
ing it together with spectra collected on other days. We 
are planning to identify and study the optimal procedure 
to subtract the ionospheric effects and systematic errors 
introduced by such a procedure. 

5.2. Power spectrum of ionospheric fluctuations 

Based on TEC measurements. IPatta et all (j20I5[) con¬ 
clude that the precision required for the detection of the 
global EoR (^mK) cannot be achieved by a ground- 
based instrument because of the flicker noise nature of 
the ionospheric fluctuations. The power of flicker noise 
is proportional to I//“, where the exponents a observed 
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in the real physical processes are in the range 0 < a < 2 

(| Pressing^ lMilottK002t lJohnsonlll925t lKirchneill2005[) . 

Therefore, the power at low frequencies can be very high, 
but not infinite as the mathematical formula suggests. 
Thus, the flicker noise does not have a well defined-mean 
value, whic h moves aw ay from its initial value as time 
progresses (lPresslll978n . 

The flicker noise of ionospheric properties would trans¬ 
late into the flicker noise of the observed sky-averaged 
spectrum preventing even a perfect ground based instru¬ 
ment from integrating down to a mK precision. However, 
it is expected that for real physical processes the flicker 
noise power law spectrum has a turn over point at some 
low frequency (/&), breaking down to white noise, and 
therefore the integral of the power does not diverge at the 
lowest frequencies. Many phy sical processe s do exhibit 
such a break in their spectra (IMilottil I2002H . In such a 
case averaging data over long periods, much longer than 
the corresponding period of the turn over point (I//?,), 
still leads to reduction in the standard error of the mean. 

Therefore, we decided to investigate the spectral char¬ 
acteristics of variations observed in the sky-averaged sig¬ 
nal (as in Fig. [S]), which we attribute to ionospheric 
variability (Sec. 13.71) . In order to extend our analysis 
to the lowest possible frequencies, we used the whole 
three months of data (~100days) at Ri35sec time res¬ 
olution, which allowed us to sample frequencies from 
Ri 0.2 X 10“®Hz (period ~49days) to r; 0.01Hz (period 
R;70sec). We only looked at nighttime data as it is the 
most important from the global EoR perspective, and 
is not contaminated by daily solar activity. The night¬ 
time power spectrum was averaged in the frequency band 
80 — 85 MHz and its LST dependence was fitted with a 
Fourier series (eq. and Fig. EOl) . Then, the Fourier 
series parameterization was subtracted from the original 
time series and the residuals are shown in Figure E51 

The resulting time series of the fluctuations has gaps 
(due to daytime, RFI excision etc.). Therefore, in or¬ 
der to determine the power spectrum of the ob served 
fluctuations we use d the Lomb-Scargle algorithm (jLombI 
Il970t [^arglelll98^ . which is suitable for unevenly sam¬ 
pled data. The power spectrum calculated by the Lomb- 
Scargle algorithm with the ~ 1// dependence is shown 
in Figure EH] Its frequency dependence at frequencies 
> 10“^ Hz (corresponding to period <1 day) can, indeed, 
be very well characterized by a flicker noise. However, 
the power law dependence breaks at the lowest frequen¬ 
cies (< 10“® Hz), suggesting a low frequency cut-off and 
therefore absence of high power contained in very long- 
period fluctuations in the data. The low frequency cut¬ 
off may explain the results of our statistical test (Sec. [5]), 
suggesting that the standard error decreases as 
even at very long integration times, as expected in the 
absence of flicker noise, and therefore makes averaging of 
multi-day (or even multi-month) data samples to ^mK 
precision possible. Analysis of even larger data sam¬ 
ples (longer time series) will enable us to verify the ob¬ 
served turn over point in the power spectrum at fre¬ 
quency fb ~ 10“® Hz. 

6. CONCLUSIONS 

Achieving the calibration precision required for the 
global EoR experiment is certainly a very challeng¬ 
ing task from the instrumental and engineering points 



Figure 25. Nighttime fluctuations of sky temperature averaged 
in the 80 — 85 MHz frequency band (residuals of Fourier series fit 
to the observed time series). 



Figure 26. The power spe ctrum of the sky temperature fluc¬ 
tuations shown in Figure (25] The frequency depenedence of the 
power spectrum exhibits flicker noise nature (~ 1//) of the varia¬ 
tions at frequencies > 10“® Hz (corresponding to periods <lday). 
However, the power law shows the evidence of a break at lower 
frequencies, which makes integration to ~mK precisions possible 
for ground-based instruments. 


of view. Moreover, besides instrumental effects there 
are propagation effects in the ionosphere, which signif- 
icantly affect the freque ncy structure of the foreground 
l|Vedantham et aP 120141 ) and therefore complicate the 
data analysis. The most relevant processes are iono¬ 
spheric absorption, emission and refraction which, based 
on variability studies of ionospheric properties (such as 
TEC), could potentia lly impede detect ion of the global 
EoR from the ground (iDatta et al.ll2015D . Therefore, it is 
important to identify the ionospheric effects in the sky- 
averaged radio spectrum at the frequencies of interest 
(50-200 MHz) and evaluate their significance for this type 
of experiment. 

In this paper, we presented analysis of three months of 
data collected between Oct 2014 and Jan 2015 with an 
upgraded BIGHORNS system equipped with a conical 
log spiral antenna and installed at the MRO. We per¬ 
formed our data analysis separately on the three indi¬ 
vidual months (Oct/Nov 2014, Dec 2014 and Jan 2015). 
We applied almost the same data reduction, RFI exci¬ 
sion and calibration procedures as described in S15 and 
extended it with a variability criterion to excise data af¬ 
fected by signihcant variability (due to for example solar 
activity or lightning). 

Our analysis was based on the assumption that the dif¬ 
ferences in spectra (averaged in 1 h LST bins) collected 
at the same LST (particularly during the nighttime) and 
similar environmental conditions (i.e. ambient temper¬ 
ature) should represent variations due to different iono¬ 
spheric conditions. Based on the sky model integrated 
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with the simulated antenna pattern, we estimated the 
expected magnitude of the ionospheric effects and com¬ 
pared the magnitude of absorption/emission against re¬ 
fraction. The expected magnitude of these effects at fre¬ 
quencies above 80 MHz is ^ 1 — 100 K. These estimations 
also led us to conclude that during the nighttime the 
contribution of ionospheric absorption and emission is 
expected to be significantly (about order of magnitude) 
larger than the contribution of refraction. The observed 
differences of the 1 h LST spectra were fitted with a sum 
of power laws Ai'~^ -|- which initially confirmed 

their ionospheric origin. 

Based on the approximate linear correlation of the fit¬ 
ted parameters A and B, which is expected if ionospheric 
absorption and emission dominate over refraction, we de¬ 
rived the electron temperature of the ionosphere as a 
function of LST and local time. The measured electron 
temperatures agree reasonably well with predictions of 
the IRI-2012 model (especially during the nighttime). 
Based on the QDC derived from 3 months of data, we de¬ 
rived the magnitude of relative absorption (optical depth 
T RiO.Ol) and its variability {6t RiO.005) as a function of 
local time, and these values are also in agreement with 
other measurements (keeping in mind they can be subject 
to differences due to geographical location, time of the 
year, or phase of the solar cycle). The analysis of a larger 
data sample should refine these results. The above re¬ 
sults lead us to conclude that we do observe ionospheric 
effects in our data and that their magnitude agrees with 
expectations. The linear correlation between the fitted 
parameters A and B indicates that, particularly during 
the nighttime, absorption and emission are indeed the 
dominant effects. 

In order to verify whether sky-averaged spectra can be 
integrated down to the required precision in the pres¬ 
ence of ionospheric effects, we have performed a statis¬ 
tical test on the three months of data. The standard 
deviation calculated from an increasing number of sam¬ 
ple spectra converges to an approximately constant value 
after the number of sample spectra is sufficiently large. 
This result is very encouraging for ground-based global 
EoR experiments because it means that the statistical 
noise can be suppressed by a factor l/^/Ns, where Ng is 
the number of averaged spectra (or alternatively Xjyjr, 
where r is the integration time). In order to verify our 
understanding of the test results we applied the same 
test to model data (without any extra variability due to 
ionosphere or instrumental effects). The standard devia¬ 
tion of the model data converges to a value slightly lower 
than observed in the the real data. The difference corre¬ 
sponds to the standard deviation of intra-day variability 
calculated from the differences of 1 h LST spectra, which 
we attributed to the ionospheric effects. Therefore, we 
conclude that the observed standard deviation can be 
understood in terms of squared sums of standard devia¬ 
tions of model data (statistical noise due to the galactic 
foreground without any ionospheric effects) and observed 
intra-day variations (attributed to ionospheric effects). 
Hence, the results of the test suggest that the spectra 
can be averaged down to the precision required by global 
EoR studies even with ground-based instruments. 

This result is in contradiction to re sults (based on TEC 
measurements and simulations) by iDatta et al.l (|2015f) 
suggesting that flicker noise (power ^ l//“) introduced 


by ionospheric variability cannot be averaged down to 
the required ~mK precision. Therefore, in order to un¬ 
derstand the apparent discrepancy between these results, 
we determined the power spectrum of the nighttime fluc¬ 
tuations observed in our data. The power spectrum 
indeed exhibits flicker noise dependence at frequencies 
> 10“^ Hz, but seems to have a turn over point (not 
uncommon in other physical processes) at a frequency 
ft! 10“® Hz (period ftl day) and seems to be constant 
(white noise) at lower frequencies. Hence, the power 
present in the lowest frequencies is finite and therefore 
long timescale fluctuations should not impede long inte¬ 
grations to suppress the standard error of the mean. We 
conclude that both results of the statistical (standard de¬ 
viation) test and power spectrum of the observed varia¬ 
tions suggest that the ionospheric stochastic error can be 
suppressed by long (multi-day or even multi-month) inte¬ 
grations. Therefore, the ionospheric absorption, emission 
and refraction and their variations are not fundamental 
impediments to integrate down the cosmic noise to the 
required ^mK precision. 

Nevertheless, because the observed ionospheric effects 
are (as expected) 1-3 orders of magnitude larger than the 
global EoR signal, they have to be taken into account 
in the data analysis if ^inK precision is to be achieved 
by a ground-based instrument in a reasonable time pe¬ 
riod. In the near future, we are planning to develop an 
appropriate procedure and study any systematic errors 
introduced. 
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